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Infinite ergodic theory is a branch of mathematics 
which deals with dynamical systems whose invariant den¬ 
sity is infinite, namely non-normalized densities describe 
statistical properties of the system. This theory was ap¬ 
plied to simple deterministic transformation rules like the 
Pomeau-Manneville map, a well known model for inter- 
mittency, where an infinite invariant density, namely a 
non-normalizable state describes long time aspects of the 
dynamics. In this context mathematicians consider two 
classes of observables, integrable and non-integrable with 
respect to the infinite density. Only recently has this con¬ 
cept gained attention in Physics (see Refs, in i)- 

Indeed at a first glance infinite densities (or more spe¬ 
cialized infinite covariant densities) might seem strange, 
since a system where the number of particles is conserved, 
must obey rigorous normalization condition for all times 
as was pointed out in the comment 0 and well known 
from elementary courses. In our work we investigated 
the Levy walk model, a widely applicable norm conserv¬ 
ing model of super-diffusion, showing that it is described 
by an infinite covariant density. Below we show that the 
comment recently published on cond-mat, is based on 
misinterpretations of the concepts and results presented 
in our work. 

To begin with, the author of the comment writes 
that we claimed that particle distribution or equivalently 
probability densities can become non-normalized in the 
case of anomalous Levy walk. This is a false accusation. 
The probability density function of the particle’s posi¬ 
tion in space, is perfectly normalized, at all times. In our 
work we used the concept of infinite covariant densities, 
i.e. non-normalized densities. However, these densities 
are not probability densities, as is well known, and clear 
from our discussion (see some details below). It seems 
that the author of the comment has assigned to the infi¬ 
nite density a meaning of a probability density, which is 
wrong. 

To further discuss the comment we first present briefly 
our main results. This is required since the author of the 


comment has not presented a full picture of our results. 
We investigated the Levy walk model and showed that 
in the super diffusive phase, the center part of the distri¬ 
bution Pcen{x,t) Eqs. (14,15) [H is described by a Levy 
stable law. This distribution is perfectly normalized, and 
a standard tool in the field. In addition we defined the 
infinite covariant density Eq. (9) 

lim t°'P{x,t) = Icd{x/t). 
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Since P{x,t) is the density, which as mentioned is per¬ 
fectly normalized P{x,t)A.x = 1 and since we have 
1 < a < 2 it is not surprising that the infinite density 
is not normalized, since the spatial integration over the 
left hand side of the equation gives infinity since ^ oo 
(this conclusion is reached by author of the comment, 
a trivial insight which is clear from Q). It is also very 
clear, even without gaining any insight on the subject, 
that the infinite covariant density is not a probability 
density function. The goal of our paper was to show in 
what sense does the concept of infinite density describe 
the statistical properties of the Levy walk model, and 
to obtain analytical expressions for this little understood 
function. First insight was that from data, e.g. numerical 
simulations of the model or in principle experiments, one 
may construct a histogram and then plot t°‘P{x, t) versus 
x/t and then observe that in the limit of large time t this 
scaled histogram approaches the analytical expressions 
for the infinite density given explicitly in [I|. This is in 
principle easy to check numerically, and it is only a pity 
that the author of the comment did not find time to do so 
(i.e., repeat simulations presented in our work i). Sec¬ 
ondly the infinite density describes high order moments, 
(|a;|‘^) with q > a. In this sense high order moments are 
integrable with respect to the infinite density, and their 
asymptotic values are obtained from this non-normalized 
density. In contrast, moments with q < a, including the 
normalization, q = 0, are non-integrable with respect to 
the infinite density, and hence they are computed with 
respect to the Levy distribution Pcenix,t). This in turn 
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is related to strong anomalous diffusion, a behavior found 
in many systems. 

We therefore found it very disturbing to read that ac¬ 
cording to the author of the comment, Eq. (8) cannot 
be applied to the whole range of x variation. As men¬ 
tioned, Eq. (14,15) in [I| explicitly give the behavior of 
the center part of the packet of particles Pcen{x,t), in 
terms of symmetric Levy distributions. Thus the center 
part of the packet is described by a Levy distribution, so 
clearly the infinite density is not applicable in the whole 
range of x. Eq. (18) gives (|a;|°) = 1, namely the normal¬ 
ization condition holds as it should. Hence the author’s 
suggestion that we claimed that probability densities can 
become non-normalizable in the ease of anomalous diffu¬ 
sion is detached from the reality of our Letter. Similarly, 
the author of the comment also writes on our Eq. (13), 
which is an equation for the moments of the process, that 
it is generally wrong, e.g. for q < a. However any one 
reading our paper sees one line before Eq. (13) that it 
is valid under the explicitly stated condition that q> a. 
Our Eq. (18) gives the solution for q < a. Unfortunately 
we see that the author of the comment did not present 
our results decently, since he does not provide the condi¬ 
tions under which our formulas work, these being clearly 
stated in [l[. 

Let us turn back to the general philosophy of the com¬ 
ment since it presents a matter of opinion, namely that 
the infinite density cannot reflect physical reality. First 
as mentioned if we scale numerical or experimental data 
as t°‘P{x, t) and plot it versus x/t the plot will approach 
the infinite density with its characteristic non integrable 
pole for small argument of T = x/t (see figures in Q)- 
In that sense the infinite density reflects physical reality. 
The claim made by the author of the comment is that 
any descent experiment, either real or numerical, done at 
finite time t will yield Jcd(u, t) which is perfectly normal¬ 
ized, and not Icd(v) and similarly Importantly, stochas¬ 
tic numerics can be done only at finite t. The general 
claim that experiments are done on finite time and hence 
asymptotic results have no value is a philosophy which 
is long abandoned. At the starting days of diffusion the¬ 
ory, one could wonder if the diffusion equation is correct. 
The solution of the diffusion equation, predicts that af¬ 


ter a fraction of a second a particle which started in New 
York city, can have a finite probability (though small) 
to be in Tokyo. Does this imply that we should throw 
away the diffusion equation, or for that aim the Gaus¬ 
sian central limit theorem, which is also valid only in a 
limit? Clearly asymptotic laws should apply within their 
limitations, and the same is true for the infinite density. 

In the case of the Levy walk and its infinite density, 
which is of-course an asymptotic result, the situation is 
even sharper. Physics actually begs for the infinite den¬ 
sity concept. The center part of the packet is described 
by the Levy stable law, as mentioned. It therefore has 
the awkward property, at first glance only, that {x'^) = oo 
at a finite asymptotic time, since the second moment of 
a Levy distribution is infinite as is well known. Thus, 
the Levy central limit theorem predicts an unphysical 
behavior that particles can travel with an infinite speed. 
Should we reject this limit theorem, because at any fi¬ 
nite time we cannot attain this divergence? In fact the 
infinite density concept cures this unphysical behavior. 
The second moment {x'^) is finite if one uses the infinite 
density concept as we have shown in [I| . We see that the 
Levy walk process has two scaling solutions. One is the 
familiar Levy density which gives a finite normalization 
and a diverging mean square displacement, the other is 
the infinite density which is not normalized but provides 
the second moment correctly. The second moment is eas¬ 
ily considered the most frequently measured moment in 
diffusion processes. Both these densities are essential for 
the correct statistical description of the process, both can 
be measured by proper scaling of data, both are strictly 
valid in the long time limit, and the practitioner in the 
field will no doubt be able to comprehend their domain 
of validity and avoid abuse. 

In summary, the authors cannot be held responsible 
if others, intentionally or unintentionally, misrepresent 
our correct results by not observing the conditions under 
which those are obtained. 
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